Unified dark energy thermodynamics: varying w and the — 1-crossing 
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We investigate, in a unified and general way, the thermodynamic properties of dark energy with 
an arbitrary, varying equation-of-state parameter w(a). We find that all quantities are well defined 
and regular for every w(a), including at the —1-crossing, with the temperature being negative in the 
phantom regime (w(a) < —1) and positive in the quintessence one (w(a) > —1). The density and 
entropy are always positive while the chemical potential can be arbitrary. At the —1-crossing, both 
temperature and chemical potential are zero. The temperature negativity can only be interpreted 
in the quantum framework. The regular behavior of all quantities at the —1-crossing, leads to the 
conclusion that such a crossing does not correspond to a phase transition, but rather to a smooth 
cross-over. 
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I. INTRODUCTION 

Recent cosmological observations support that the uni- 
verse is experiencing an accelerated expansion, and that 
the transition to the accelerated phase has been realized 
in the recent cosmological past [1]. In order to explain 
this remarkable behavior, and despite the intuition that 
this can be achieved only through a fundamental theory 
of nature, physicists can still propose some paradigms for 
its description, such are theories of modified gravity , 
or "field" models of dark energy. The field models that 
have been discussed widely in the literature consider a 
canonical scalar field (quintessence) [H , a phantom field, 
that is a scalar field with a negative sign of the kinetic 
term or the combination of quintessence and phan- 
tom in a unified model named quintom [f| . The common 
feature of all the field models of dark energy, is that the 
equation-of-state parameter w of the dark energy "fluid" , 
is induced by the field evolution and thus it acquires a 
dynamical nature (w — > w(a) with a the scale factor). 

On the other hand, there have been large efforts in 
order to study the thermodynamic properties of dark en- 
ergy fluids with a constant w = Wq 5^ — 1. For the case 
of quintessence scenarios, where w = wo > — 1, one can 
use the results of conventional perfect fluid thermody- 
namics. However, in the case of phantom fields, where 
w — wo < —1, a thermodynamic investigation must be 
performed from first lines. Some authors assume a zero 
chemical potential fj, and argue that the temperature of 
a phantom fluid must be negative, with the density and 
the entropy positive @, 0], and the negative nature of 
the temperature is related to the quantum properties of 
phantom fields (although in some initial works on /i = 
the entropy was found to be negative and thus the phan- 
tom phase meaningless [8]). On the other hand, some au- 
thors consider the non-zero chemical potential case, and 
argue that phantom fluids must have /i < 0, with tem- 
perature, entropy and density being positive, with the 
phantom particles possessing a bosonic, massless nature 



In our opinion, the central point of dark energy models 
is not so much that w can take a constant value different 
from —1, but rather that w(a) is varying (indeed since 
w(a) is a ratio of two varying quantities, being a con- 
stant requires rather special constructions). In addition, 
the crossing of the phantom-divide w = — 1 from above 
to below, is not only possible but it could indeed be the 
basic requirement for a successful description of observa- 
tions. Therefore, the question of investigating in a unified 
way the thermodynamics of dark energy fluids with an 
arbitrarily varying w(a), becomes crucial. 



However, a comment must be made here, concerning 
the applicability of thermodynamics to time-dependent 
systems such as the universe, where gravitational degrees 
of freedom are present. Although in principle statisti- 
cal mechanics assumes equilibrium and ergodicity, which 
could be violated in time-dependent gravitational back- 
grounds, and despite the fact that a thermodynamic and 
statistical description of the dynamics of gravitational 
degrees of freedom does not fully exist yet, we believe 
that such an approach can be enlightening. This is sim- 
ilar to the thermodynamical investigation of the prop- 
erties of black holes or de Sitter spaces, that is purely 
gravitational systems evolving in time. In addition, it is 
fairly known that black holes can be regarded as the time- 
reversed analog of the Big Bang universe, despite start- 
ing from static spacetime metrics. Having these points 
in mind we proceed to the thermodynamic analysis of a 
general dark-energy fluid. The paper can be outlined as 
follows: In Sec. [IT] we present the thermodynamic prop- 
erties of a general-w;(a) dark-energy fluid, while Sec. IIIII 
is dedicated to discussion, concluding remarks and a brief 
summary of the results. 
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II. GENERAL w(a) THERMODYNAMICS 

We consider dark-energy fluids, with energy density 
and pressure respectively p(a) and p{a), and with an ar- 
bitrary and general equation of state of the form: 



p(a) = w(a)p(a). 



(1) 



In these expressions a is the scale factor of a homoge- 
nous and isotropic Friedmann-Robertson- Walker uni- 
verse, with metric: 



ds 2 = dt 2 -a 2 (t) 



dr 2 



1 — kr 2 



r 2 d6 2 



r 2 sin 2 i 



(2) 

where k = 0, ±1 is the curvature parameter and t is 
the comoving time. Since we desire our analysis to be 
general, we do not consider specific cosmological mod- 
els, with explicitly extracted Friedmann equations, but 
we just consider the resulting 10(a). Thus, the present 
investigation is valid for every FRW cosmological model. 
Finally, note that we could equivalently express w(a) as 
w(t) or w(z) (with z the redshift), since such a transfor- 
mation is straightforward. 

As usual, the equilibrium thermodynamic states of a 
relativistic perfect fluid are characterized by an energy- 
momentum tensor T^ v : 



= (p + p)u^u u - pg" 



(3) 



where (p ~ {0, 1,2,3}) is the 4-velocity in the met- 
ric g^" . The particle and entropy currents, and 
respectively, are defined as: 



AT" 



(4) 



with n and s the densities of particle number and entropy. 
Denoting the covariant derivative by (;), the equations 
of motion are given by T^ v \ v — 0, iV% and S%= 0, 
which in the case of the FRW geometry write: 



p + 3[l + w(a)}p- 



a 

a 
>,— 
a 

a 






(5) 





(6) 


0. 


(7) 



3s- 



with the dot denoting the comoving time derivative. In 
the case of a general w(a), the solution for p(a) reads: 



p(a) = p exp 



10 3[w(a) + 1] 



da 



p(a) = po 



„3[to +l] 

"0 

2 3[™(a) + l] 



exp 



da w'(a) In a 



,(8) 



where the prime denotes the a-derivative, and ao and 
Wo = w(a ) are the present values of the scale factor 
and of the dark energy equation-of-state parameter re- 
spectively. In this relation we preferred to use u>'(a), in 



order for the extraction of the w = const results (i.e. 
w'(a) — » 0) to be straightforward. Similarly, the solu- 
tions for n(a) and s(a) are simply: 



n(a) = n ^) , s(a) = s ^— j 



(9) 



with no and So the values of the corresponding quanti- 
ties at present (from now on the index stands for the 
present value of a quantity). Finally, we stress that w(a) 
is the value of the equation-of-state parameter when the 
scale factor is a, while wq is its present value. Obviously, 
on can face all the evolution combinations, where w(a) 
can cross —1 at a specific scale factor or not, and/or wq 
being —1 or not, or the case w(a) — wq — const, with 
the constant being not equal or equal to —1 (the last 
case is the cosmological constant universe). Therefore, 
in the following one must examine the regularity of all 
expressions when w(a) and/or wq is equal to —1. 

Let us now refer to thermodynamics. Considering 
(n, T) as the 2D thermodynamical base, we use the iden- 

, which in 

the case of a perfect fluid writes: 



tity (Gibbs law): T | = p + p 



T 



dw(a)p(a) 



= [w(a) + l]p(a)—n(a) 



dp{a) 
dn 



(10) 



dT 

On the other hand, we can express p as p = 
n (f^) + T (jf) • Combining this relation with (TO] 

and (5|) we acquire 3| T(a)w(a) + T(a 
-3^T(a)p(a)w'(a), or using (SJ): 





dp(a) 




da 



-T{a)w(a)+T'(a) 
a 



[w(a) + l] = T(a)w'(a). (11) 



Note that this general differential equation leads straight- 
forwardly to the general physical implication that T(a) 
must be zero when w(a) crosses —1. Furthermore, note 
that in the case where w(a) = —1 = const, at all times, 
i.e. of a pure cosmological constant, equation pip is triv- 
ially satisfied, not allowing for the determination of T(a), 
and this has led some authors to argue that the temper- 
ature of the vacuum state is ill-defined [7|. However, we 
see that considering a general w(a), i.e with w'(a) =/= 0, 
we are led to T = for the vacuum state, which is a 
self-consistent result. 

If w(a) 7^—1, the general solution of (| 1 1 [) is: 



T(a) = To 



w(a) + 1 



Wo + 1 



exp 



-3 



daw' (a) In a 



(12) 

with To the present temperature value. Solution (T2|) 
holds for w(a) < —1 or w(a) > — 1. For w(a) = —1 the 
only information that we can obtain straightaway from 
(jlip is that T(a) must be zero. However, this requirement 
is indeed satisfied by (112|) , and furthermore we can see 
that: 



lim T(a) 



0. 



(13) 
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Thus, we can safely consider that (| 1 2[) is the general solu- 
tion for every w(a), including the — 1-crossing. We men- 
tion that if wo — > — 1, then as required by (jlljl To — ► 
too, and T(a) remains regular. 

Expression (|12[) has an interesting and general physical 
implication, namely that T(a) and w(a) + 1 (and thus To 
and wo + 1) have always the same sign. Therefore, in sum- 
mary, the temperature is always negative for w(a) < — 1, 
it is zero for the cosmological constant bound, and it is 
always positive for w(a) > —1. Later on, we will discuss 
these points in detail. Finally note that we could denote 
the regular ratio To/(wo + 1) by a positive constant Co, 
but we prefer keeping it in order for the various quantities 
at present to be straightforwardly obtained. 

A useful relation can be obtained using ([5]) and (p"2")l . 
namely: 



P(a) = Po 



T{a) [wp + X\ 
T [w(a) + l] 



(t) • 



(14) 



which is regular at w(a) = —1 (and at wq = —1), as 
can be seen from (TT2"1) . Equivalently, eliminating a 3 be- 
tween (|T2|) and (TH]) we acquire the generalized Stefan- 
Boltzmann law: 



P{T) = Po 



T(q) [wq + 1] ) =ra 
T [tu(o) + l] 



Hal 



• exp 



«;(a) 



daw' (a) In a 



(15) 



Note that in the case w(a) = 0, i.e when the dark energy 
fluid is a dust, one can only use relation (I14|) . since the 



elimination of a between (|12|) and ([14]) is not possible. 

Let us now return to thermodynamics, using the gen- 
eral identity Ts — p + p — pn, with p the chemical po- 
tential of the dark energy fluid, considered non-zero in 
general. This relation can be used to determine p. In 
particular, written as 



p(a) 



[w(a) + l]p(a) T(a)s{a) 



n(a) 



n(a) 



(16) 



and using ((8|) and (fT2|) it leads to 



p(a) = fi 



w(a) + 1 



w + 1 



a 



~.3w(a) 



exp 



r a o 

-3 / daw' (a) In a 

J a 



where we have defined 



Mo 



_ Pq(wq + 1) - TqSq 



n 



(17) 



(18) 



the chemical potential value at present. Expressions (fT7|) 
and (fT8|) are regular at wo = — 1 (since in this case To = 
too), and furthermore at w(a) = — 1 (TTT)) leads to /j,(a) = 
0. Note however that in general the sign of po, and thus 
of /x(a), can be arbitrary. Therefore, we conclude that the 
chemical potential for the cosmological constant bound 
is zero, but it is arbitrary at the two sides of that bound, 



depending on the specific cosmological model. Later on, 
we will discuss these points in detail. Finally, (fTB|) with 
the use of lfl4|) . leads to the simple relation: 



P(T) = Po 



T(a) 



To 



(19) 



Relation (|T6|) . with s(a) — S(a)/V(a) 7 can be used 
to calculate the total entropy 5* of the universe, with 



physical volume V(a) 
([!!))) ■ we acquire: 



Indeed, using ©, (fl~5j) and 



S(T) = s V(a) 



T(c 



[w 



T [w{a) + l] 



■ exp 



w(a) 



3[w(a)—w ] 

a n l ° (a) 



da w'(a) In a 



(20) 



which is regular and non-zero at w(a) 
easily seen, using (TJJ}, that 



T. It can be 



S(a) = s V(a) f ^ 
V a 



sqVq = s(a)V(a) 



(21) 



with Vo — &q the current physical volume of the universe. 
Expression (|21[) is just the statement of the conservation 
of entropy in the whole universe, and provides a self- 
consistency test for our calculations. Finally, note that 
in the case w(a) = 0, i.e when the dark energy fluid is a 
dust, relations (|20|) and (|15|) lead to a trivial result, and 
one can use only relation (f2"Tj) . 

In the aforementioned investigation we desired to re- 
main general, and we did not use any ansatz for uu(a). 
However, phenomenologically, one can consider vari- 
ous w(z)-parametrizations [ll| at will, in order to ex- 
tract quantitative predictions for the quantities at hand. 
Lastly, taking the limit w'(a) — > (or equivalently just 
setting w'(a) = since the terms that need caution have 
been written as separated pre- factors) , we extract the 
corresponding relations for the w(a) = wq = const case: 

PW = Po [—) 

/a \ 3w ° 
P(a) = Po {—) 



P(T) = Po 
P(T) = Po 



a 

T(a) 
To 

T(a) 
To 



-i+i 

"0 



S(T)=s V(a){^ 



(22) 



We mention that these relations are valid also for the 
simple cosmological constant (u>o = — 1 and Tq = n a — 
0), giving p = po = const., T = p = and S — sqVo = 
s(a)V(a) = const.. 
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III. DISCUSSION 

Let us discuss about the physical context of the ob- 
tained results. First of all, we see that all quantities 
are regular and well-defined, for all values of w(a). In 
addition, the density and entropy are always positive, 
consistently with the basic requirements of classical and 
quantum field theory. This behavior clarifies some am- 
biguities about the phantom nature, since one does not 
need integrability conditions Q or special constructions 
[1, M, El in order to make this "phase" possible (which 
remarkably might indeed be the current phase of the uni- 
verse). 

An important physical consequence is that the tem- 
perature of a dark energy fluid with w(a) < —1 is nega- 
tive, independently of the value of the chemical poten- 
tial which can be arbitrary. This is in contrast with 
the works that required necessarily a negative /x and a 
positive T for the constant-u> phantom fluids 0, IToj |. 
The misleading point in these considerations were that 
since 10(a) = wq = const, the authors instead of equa- 
tion dTTJ) used 3T(a)w(a) + aT'(a) = 0, which leads to 
qualitative different results. The reason for this behav- 
ior is that we face a singular perturbation problem [T3 |. 
i.e the results for w'(a) — > are more general and do 
not coincide in the whole variable-range with those for 
w(a) = wo = const.. Furthermore, note that if one con- 
siders w(a) = wq = const, then he cannot educe any 
result at all for the temperature sign, since this will coin- 
cide with the sign of To, which is arbitrary in principle. 
Only allowing for a — 1-crossing (i.e for a general w(a)) 
we can safely conclude that T(a) and w(a) + 1 (and thus 
To and wq + 1) have the same sign, as can be seen from 

Another significant physical implication, is that the 
vacuum, either as a permanent state (w(a) = wq = — 1 
at all times) or as an instantaneous state at the time of 
the crossing of a varying- w (a), has zero temperature and 
zero chemical potential. This was expected, since we can 
consider it to correspond to the zero mode of the various 
involved scalar fields. On the contrary, a vacuum with 
non-zero T and /x would require a non-trivial explanation. 
We mention that starting with w(a) = wo — —1, one 
cannot find a relation for the temperature (and thus for 
the chemical potential too), and therefore he concludes 
that the temperature and chemical potential for this state 
is ambiguous or ill-defined 7] , or he is led to wrong results 
@) EH ■ The correct approach is to start with a general 
w(a) and then examine the limit w(a) —* — 1. 

The regular behavior of all quantities at the phantom- 
divide crossing, leads to an additional interesting result. 
In particular, we conclude that such a crossing does not 



correspond to a phase transition, but rather to a cross- 
over. An additional argument for this statement arises 
from the regular behavior of the "specific heat" , defined 
as cx ^ or cx T^y, at the crossing. The —1-crossing, 
dynamically is just a pass to the super-accelerated evolu- 
tion, while thermodynamically is just a smooth pass from 
T > to T < 0, without a phase transition. This behav- 
ior was observationally required, since if a radical cos- 
mological phase transition had taken place in the recent 
cosmological past, it would have left observable imprints. 

The temperature negativity can only be interpreted in 
the quantum framework. Although the discussion about 
the construction of quantum field theory of phantoms is 
still open in the literature (see for example [l3[ about 
causality and stability problems and the possible spon- 
taneous breakdown of the vacuum into phantoms and 
conventional particles), there have been serious attempts 
in overcoming these difficulties and construct a phantom 
theory consistent with the basic requirements of quantum 
field theory [3]. Thus, as it was analyzed in detail in [y], 
in such a quantum consideration of the phantom fluid 
all novel phenomena (stimulated absorption of phantom 
energy, generalized Wien and Planck radiation laws) can 
be naturally positioned. 

One of these novel phenomena, that can have inter- 
esting cosmological implications, is the accretion of dark 
energy onto black holes. In particular, in [TEl ] it was 
shown that the solution for a stationary, spherically sym- 
metric, accretion of a cosmological perfect fluid onto 
a Schwarzschild black hole with mass M, is given by 
M = 4irAM 2 (p +p), with A a positive constant, which 
in the case of a general w(a) reads: 

M = 4ttAM 2 [1 + w(a)]p(a). (23) 

Although a specific solution of this equation requires a 
constant w(a) (see @), in general we observe that M > 
in the quintessence regime, M < in the phantom one, 
while M = at the —1-crossing, consistently with the 
generalized second law of thermodynamics (the Hawking 
radiation is neglected in this consideration). Relation 
(I2U1) is independent of the phantom chemical potential 
(contrary to the requirement of [16]) and the change in 
the sign of M is consistent with the change in the T-sign 
@. Although smooth at the —1-crossing, such a change 
in the variation rate of a black hole mass, and thus of 
the total mass in black holes in the universe, could leave 
observable imprints in the celestial orbits. 
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